o 



O 

o 
> 

•i-H 

X 



Simulating pasta phases by molecular dynamics and cold atoms — 
Formation in supernovae and superfluid neutrons in neutron stars 
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In dense stars such as collapsing cores of supernovae and neutron stars, nuclear "pasta" 
such as rod- like and slab- like nuclei are speculated to exist. However, whether or not they 
are actually formed in supernova cores is still unclear. Here we solve this problem by demon- 
strating that a lattice of rod-like nuclei is formed from a bcc lattice by compression. We also 
find that the formation process is triggered by an attractive force between nearest neigh- 
bor nuclei, which starts to act when their density profile overlaps, rather than the fission 
instability. We also discuss the connection between pasta phases in neutron star crusts and 
ultracold Fermi gases. 



£^j . §1. Introduction 

Collapse driven supernova explosion, an explosion in the death of a massive star, 
is one of the most dramatic phenomena in the universe and has been a long-standing 
mystery in astrophysics. X) One of the key ingredients to understand the mechanism 
of the supernova explosion is the study of matter in the core of supernovae. 2) 

State of matter in supernova cores also undergoes dramatic changes in the pro- 
cess of the collapse. Matter experiences an adiabatic compression, in which the 
density in the central region of the core increases from ~ 10 9 g cm -3 in the begin- 
ning of the collapse and finally reaches around the normal nuclear density ^3x 10 14 
g cm -3 (corresponding to the number density of nucleons po = 0.165 fm -3 ) just 
before the bounce. It is predicted that, in the final stage of the collapse, nuclei 
are rod-like or slab-like rather than (roughly) spherical in the central region of the 
coreP™ These non-spherical nuclei are collectively called nuclear "pasta". 

Recently, nuclear pasta attracts much attention of astrophysicists and nuclear 
physicists [see, e.g., Refs. ED®] since it has been pointed out that the nuclear pasta 
can occupy 10 - 20 % of the mass of supernova cores^ and thus might have influences 
on the dynamics of supernova explosions.^"® However, a fundamental problem of 
whether or not and how pasta phases are formed during the collapse is still unclear. 
In the present workp^J we succeed in numerically simulating the formation of a 
trianglar lattice of rod-like nuclei from a bcc lattice of spherical nuclei in collapsing 
cores and solve the above long-standing problem. In addition, we discover that 
the formation process is very different from a generally accepted scenario based on 
an instability with respect to nuclear fission. Our work provides a solid basis for 
understanding materials in supernova cores, which is indispensable for solving a 
long-standing problem of the mechanism of supernova explosions. 
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§2. Fission Instability? 

In a generally accepted conjecture, based on the Bohr- Wheeler condition for 
vanishing the fission barrier derived for isolated nuclei, it is predicted that the for- 
mation of the pasta phases are triggered by the fission instability with respect to the 
quadrupolar deformation of spherical nuclei P^* The essential point of this prediction 
is that, at higher densities, the effect of the Coulomb repulsion between protons in 
nuclei, which tends to make a nucleus deform, becomes dominant over the effect of 
the surface tension of the nuclei, which favors a spherical nucleus. However, Refs.[T2]) 
and [13]) have pointed out that the background electrons, which have been ignored in 
the above prediction, suppress the effect of the Coulomb repulsion between protons 
in nuclei and the fission barrier never vanishes in the relevant density region. These 
findings cast a doubt on the above conjecture based on the fission instability. 

§3. Theoretical Framework 

Since formation of the pasta phases is accompanied by dynamical and drastic 
changes of the nuclear structure, an ab-initio approach is called for. The quantum 
molecular dynamics (QMD), 14) which can properly incorporate the thermal fluctua- 
tions^ and enables us to simulate large systems with ^ 10 3 nucleons or more™ is 
a suitable approach for the present purpose. 

We use the QMD Hamiltonian of Ref. 17). In our simulations, we consider a 
system with protons, neutrons, and charge-neutralizing background electrons in a 
cubic box with periodic boundary condition. Here, we shall focus on the case of the 
proton fraction x ~ 0.39 (total number of nucleons N = 3328 with 1312 protons and 
2016 neutrons). 

We simulate the compression of the bcc phase of spherical nuclei in the collapse 
[see Ref. [T0]1 for detailed procedures]. Starting from p = 0.15/>o 5 we increase the 
density by changing the box size L (the particle positions are rescaled at the same 
time). Here the average rate of the compression is < O(10~ 6 ) po/(fm/c) yielding the 
time scale of > 10 5 fm/c to reach the typical density region of the phase with rod- 
like nuclei. This is much larger than the time scale of the change of nuclear shape 
and thus the dynamics observed in our simulation is determined by the intrinsic 
physical properties of the system. We perform adiabatic compression and isothermal 
compression at various temperatures [see Ref. [TU]) for details]. In all the cases, we 
observe the formation of rod-like nuclei; here we show a typical example in which we 
obtain a clear lattice structure of the rod-like nuclei. 

§4. Results 

Figure Q] shows the snapshots of the formation process of the pasta phase in 
adiabatic compression. Here, we start from initial condition at p — 0.15po and 
T = 0.25 MeV (t = fm/c) [Fig. [Q(a)]. At p ~ 0.243p [Fig. 11(c)], the first pair of 
two nearest-neighbor nuclei start to touch and fuse (dotted circle), and then form 
an elongated nucleus [see, e.g., Fig. 11(d)] . After multiple pairs of nuclei become 
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Fig. 1. Snapshots of the transition process from the bcc lattice of spherical nuclei to the pasta 
phase with rod-like nuclei. The red particles show protons and the green ones neutrons. In 
panels (a)-(g) and (h-1), nucleons in a limited region relevant for the two rod-like nuclei in 
the final state (h) [surrounded by the dotted lines in panel (h-2)] are shown for visibility. The 
vertices of the dashed lines in panels (a) and (d) show the equilibrium positions of nuclei in the 
bcc lattice and their positions in the direction of the line of sight are indicated by the size of 
the circles: vertices with a large circle, with a small circle, and those without a circle are in 
the first, second, and third lattice plane, respectively. The solid lines in panel (d) represent the 
direction of the two elongated nuclei: they take zigzag configuration. The box sizes are rescaled 
to be equal in the figures. This figure is taken from Ref. 10). 



such elongated nuclei, we observe zigzag structure as shown in Fig. H^d). Then these 
elongated nuclei stick together [see Figs. [Jje) and (f)], and all the nuclei fuse to 
form rod-like nuclei as shown in Fig. [U(g). Finally, we obtain a triangular lattice of 
rod-like nuclei after relaxation [Figs. [TJ h-1) and (h-2)]. 

Note that before nuclei deform to be elongated due to the fission instability, they 
stick together keeping their spherical shape [see Fig. U[c)]. Besides, in the middle of 
the transition process, pair of spherical nuclei get closer to fuse in a way such that 
the resulting elongated nuclei take a zigzag configuration and they further connect 
to form wavy rod-like nuclei. This feature is observed in all the other cases in which 
we obtain a clear lattice structure of rod-like nuclei, and the above scenario of the 



4 Gentaro Watanabe 

transition process is qualitatively the same also for those cases. It is very different 
from a generally accepted picture that all the nuclei elongate in the same direction 
along the global axis of the resulting rod-like nuclei and they join up to form straight 
rod-like nuclei P* 

Instead of such a scenario, we have shown that the pasta phases are formed in 
the following process [see Ref. [TO]) for details]. When nearest neighbor nuclei are so 
close that the tails of their density profile overlap with each other, net attractive 
interaction between these nuclei starts to act. This internuclear attraction leads to 
the spontaneous breaking of the bcc lattice and triggers the formation of the pasta 
phases. 

§5. Simulating Pasta Phases by Ultracold Fermi Gases 

In the remaining part of this article, we shall discuss the connection between the 
pasta phases and the ultracold atomic Fermi gases. Pasta phases can exist also in 
crusts of neutron stars. There, pasta nuclei are immersed in background electrons 
and a gas of dripped neutrons, which is regarded to be in a superfluid state. 

It has been pointed out that the low density neutron matter shows the same 
universal property as unitary Fermi gases provided that the interparticle separation 
r s between neutrons is much smaller than the absolute value of the scattering length 
a ~ —18.5 fm and is much larger than the effective range r e ~ 2.7 fnflH^LEj) [ see a i so 
Refs. [20]) — [22]) ] . Here, "universal" means that properties do not depend on details of 
the nuclear potential. 

In dripped neutron gas in the pasta phases, the density is p n = fcp/(37r 2 ) > 
0.05 fm -3 and thus the interparticle separation r s > 1.7 fm is comparable to r e 
(here, fcp is the Fermi wave vector of the ideal Fermi gas with the same density). 
However, Schwenk and PethidP^'^ have shown that in such a region of /cpr e <~ 1, 
corrections from the universality in the equation of state is accurately described 
within the effective range expansion and no further details of nuclear forces are not 
needed. The energy per particle E/N can be expressed in the same form as the 
universal case but with a r e -dependent factor £(/cFr e ): E/N = £,(kYr e ) (3/5)£?f 
with Ey = T?k\/2vn and m being the fermion (neutron, in the present case) mass. 
This means that the dripped neutron gas in the pasta phases can be simulated by 
cold atomic Fermi gases with a narrow Feshbach resonance. At qualitative level, we 
can also expect that unitary Fermi gases would provide a suggestive guidance for 
exploring interesting physical effects related to the dripped superfluid neutrons in 
the pasta phases. 

In Refs. [23]) and [25]) . we have studied superfluid unitary Fermi gases in a one- 
dimensional (ID) periodic potential V ex t(z) = sE^sm 2 (q^z), where s is the laser 
intensity, i?R = h 2 q^/2m is the recoil energy, q^ = 7v/d is the Bragg wave vector, 
and d is the lattice constant. This setup resembles superfluid neutrons in the pasta 
phase with slab-like nuclei, where p n ~ 0.5po — 0.08 fm -3 and d ~ 15 - 20 fm 
and thus Ey/Er = (/cf/^b) 2 ^ 40 - 70. The strength Vq = sEr of the optical 
lattice corresponds to the depth of the bottom of the conduction band of neutrons 
measured from the average potential energy outside nuclei. Since the neutron Fermi 
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Fig. 2. (Color online) Inverse compressibility ft -1 , and effective mass m* of the unitary Fermi 
gas for s = 1 (red), 2.5 (blue), and 5 (green). The solid, dashed, and dash-dotted lines show 
the results obtained by the Bogoliubov-de Gennes (BdG) equations. The dotted lines shows 
the asymptotic expressions [Eqs. (|5-1|) and (|52|) with r e = and d n £, — 0] obtained by the 
hydrodynamic theory. The s — 1 results for m* are also shown in the inset in the linear scale. 
Adapted from Fig. (1) in Ref. l2l)) . 



energy I?F,n for the density of neutrons inside nuclei is Ef jU ^35-40 MeV and 
the neutron chemical potential ji n inside nuclei is ji n ~ 10 - 15 MeV, we estimate 
Vo ~ \fJ>n — ^F,n| ~ 25 MeV. Thus, we obtain 5 ~ 25 - 50. 

In Fig. [21 we show the incompressibility ft -1 = n<9^e and the effective mass 



m 



n 



(d% 



of the unitary Fermi gas. Here e = e(n, P) is the energy density 
averaged over the unit cell as a function of the average (coarse-grained) density n 
and the quasimomentum P of the superflow in the z direction. 

The rapid reduction of ft -1 with decreasing Ey/E^ at small Ey/E^ is due to 
the formation of bosonic molecules induced by the periodic potential. Maximum of 
ft -1 and TTi* around Ey/E^ ^ 1 is the effect of the energy band gap. Since the 
tunneling rate through the barriers, which relates to m*, depends on the barrier 
height exponentially, m* increases drastically for larger s. The enhancement of m* 
has been found also in dripped neutron gas in the crust of neutron stars P^ 

Taking account of the density dependence of £ in the case of the large effective 
range, we derive k~ 1 and m* for sE^/Ey <C 1 using the equation of state, E/N — 
(3/5)£I?f, within the hydrodynamic theory. In the relevant region of /cpr e > 3, d^ is 
negligible, and keeping up to the first order of d n ^ (n^~ 1 d n ^ > < 0.1 in this regiorP^), 
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we obtain 



r I 32^ 



^(f) 2 + 3„ r ^f 1 ^(t) 



and 



m ' 32^ \E F J 



(5-1) 



(5-2) 



Since £ is a monotonically increasing function of /cp^e? m * f° r non-zero r e is smaller 
than that of unitary Fermi gases. 
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